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Increasing demands for reliable and least-cost high- 
performance aerostructures are pressing design analyses, 
materials, and manufacturing processes to new and narrowly 
experienced performance and verification technologies. This 
study assessed the adequacy of current experimental 
verification of the traditional binding ultimate safety factor 1 
which covers rare events in which no statistical design data 
exist. Because large high-performance structures are 
inherently very flexible, boundary rotations and deflections 
under externally applied loads approaching Fracture may 
distort their transmission and unknowingly accept 
submarginal structures or prematurely fracturing reliable 
ones. A technique was developed, using measured strains 
from back-to-back surface mounted gauges, to analyze, 
define, and monitor induced moments and plane forces 
through progressive material changes from total-elastic to 
total-inelastic zones within the structural element cross 
section. Deviations from specified test loads are identified by 
the consecutively changing ratios of moment- to- axial load. 

Nomenclature 

N = normal load, kips 

M = moment, kip-inches 

C = cross sectional limits, inches 
H = element thickness, inches 

w = clement width, inches 

E = elastic modulus, ksi 

n = strain hardening exponent 

K = strength coefficient, ksi 

cr = normal stress, ksi 

£ = normal strain 

Subscripts 

ty - tensile yield 

t u - tensile ultimate 

cy - compression yield 

N = normal 

M = bending 

1 = minimum measured strain 

2 = maximum measured strain 

k - zone number 
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Increasing demands for more reliable and affordable 
access to space are promoting leaner and more innovative 
structural designs that invoke more reliance on experimental 
verification of their behavior and safety . This compelling 
shift raises concerns on how well verification tests are 
implemented. Assemblies of large, high-performance 
aerostructures are inherently very flexible, and structural 
boundary rotations and deflections at externally applied loads 
approaching rupture may improperly transmit the binding 
verification loads and unknowingly reject a perfectly adequate 
design or accept a submarginal one. 

To sample this phenomenon, the slope and 
deflection were calculated at the free -end boundary 7 load on a 
hypothetical cantilevered beam, Fig. 1. Deflections and 
slopes are shown for the yield and ultimate strain limits 
calculated at the fixed end. The tangent of the free boundary 
slope 0 is a measure of the consecutively applied load 
decomposing from bending to bending-axial load ratio. 
Though the vertical scale is exaggerated, the slopes and 
displacements at the free end are proportional. 



Fig. 1. Boundary load deformation 

The ultimate load to fracture was calculated to be 
twice the yield load, and the resulting ultimate strain at the 
fixed end was an order-of-magnitude larger than the yield. At 
ultimate loading, the predicted deflection was an 18 degree 
slope resulting in over 30 percent bending-to-axial load ratio 
deviation. The adversity of this ratio to the verification 
criteria is dependent on how it feeds into and intensifies 
critically stressed regions and how it may change the failure 
mode. 

Though many codes and texts are available for 
predicting inelastic strain responses from imposed inplane 
and bending loads, literature is mute on determining test 
combined loads from measured inelastic strains. A technique 
was developed to analyze elastic-inelastic strains measured 
from back-to-back surface mounted gauges to determine and 
verify transmitted test loads with specified applied loads. 

II, Elastic-inelastic materials model 

Modeling elastic-inelastic behavior could be very 
difficult^ unless idealized into the simplest mathematical 


expressions wiLhin the physical phenomena of the material 
and its application, such as the two parameter power 
expression,-* 


a =K e" , (l) 


where “n” is the strain-hardening exponent. In the linear 
elastic region, <j<F(y. the exponent is defined as n = 1.0, 
and for cj> Fy , the strain-hardening exponent is calculated 
from uniaxial stress-strain data 


log (FJFJ 

log (£ tu / £„) ' 


( 2 ) 


The strength coefficient “K" is evaluated at the yield stress, 
which is the elastic-inelastic interface, 

K=^-=E n F^ ) . (3) 

£ty 

These properties are directly applicable to normal stresses 
and strains without interpretation through theory. 

III. Structural Modeling 

A rectangular cross section element illustrated in 
Fig 2 represents most structural components and regions as 
in beams, plates, and shells. Bending and inplane normal 
loadings are the most commonly measured components on 
this type of element using back-to-back strain gauges. They 
often may be sufficient to sample and verify the load 
transmission of more complex systems, including transverse 
shears. 


Surface mounted 
strain gauges 

1 r 



Fig. 2. Back-to-back instrumented element 


bending and axial strains are linear, they may be 
algebraically added as shown in figure 3(a). These combined 
strains are measured back-to-back at the element surfaces as 
£2 and where £2 is assumed to be greater than ej . Fig. 
3(b) illustrates the nonlinear bending stress distribution 
derived from the strain distribution using Eq. (1) and the 
shift of the bending axis to balance the moment. 



Fig. 3. Strain and stress distributions along the 
element cross section 

When £ 2 < e ty , the combined strain distributions 
over the cross section are all in one elastic zone, and when 
£ ; > E t y , the distributions are all in one inelastic zone. The 

objective is to define the elastic-inelastic zone boundaries for 
all other measured strain combinations and to calculate their 
contributions to the total normal and bending loads. 

The net strain from any miriplanc y-distance along 
the element thickness in figure 3(a) is defined by the 
proportionality 

e y =7(0.5 H +y) +£ { , (6) 

and its location is 

V =-t(^-£ 1 )-0.5 H . (7) 


The induced normal stress in Fig. 2, 


The bending strain slope is 


and strain, 



are uniformly distributed over the element cross section 
throughout the elastic and inelastic range. Because cross 
section planes are known to remain plane after elastic and 
inelastic bending, the inelastic bending strain also varies 
linearly along the thickness. However, the stress varies 
nonlinearly with Eq. (1) and the bending neutral axis is not 
expected to coincide with the cross section centroid. Since 




H 


( 8 ) 


The incremental normal load along the cross 
section thickness is the product of the induced stress and 
area, 

dN = w Gy dy - wK[e$dy . 

Substituting Eq. (6) for the strain and integrating, all zone 
normal loads may be calculated from 



(9) 


( 16 ) 
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where C a and Cb are the integration limits of a zone. A 
zone is bound along the y-axis by the surface measured 
strains, £/ and e 2f or by the material limit changes noted by 
s^and e c> . Substituting the appropriate pair of boundary 
strains into Eq. (7), 


Ca,b y ^0 


K 

1 


( 10 ) 


provides the upper and lower integration limits of each zone. 
The yield strain may be tension or compression, where £ Cy 
~ _ £ is assumed for a symmetrical material. The normal 
load across the thickness is the sum of all the zone normal 
loads 


n =2X. (id 

Bending strain along the thickness is given by e My 
= £y - and the neutral bending axis is defined by a zero 

bending strain {s My = 0 ), from which £y = £^ . 

Substituting into Eq. (7), the neutral bending axis is 

C M = | (£*-£,) -0.5 H , (12) 

where the normal strain, e N , across the thickness is 
determined by substituting Eq. (11) into Eq. (5). Using Eq. 
(12), the incremental moment about the neutral axis is 

dM - w< 7 y (y - C M ) dy = wK{eJf(y - C M ) dy . 

(13) 

Substituting Eq.s (6) and (12) into Eq. (13) and integrating, 
a zone moment about the neutral axis is calculated from 


M k =wKy 



(14) 


The moment about the cross section is the sum of all the 
zone moments. 


M = £M k . (15) 

A unit width, w = 1, is assumed for plates and shells from 
which normal loads and bending moments are defined by 
kips per inch and kip-inch per inch units, respectively. 
Using the strain distribution expression of Eq. (6), the stress 
distribuuon along each zone is given by 


<j y = K [ABS (E y )] n SGN (£ y ) 

Expressions in absolute form allow raising strains to odd 
powers. SGN ( ) is the signum function, which re- 
establishes the sign of the expression. If the function equals 
-1, the strain is negative. 

IV. Normal Load and Momen t Solutions 

As the induced normal and bending loads in figure 2 
increase, the strain distribution over the element cross 
section progresses from totally elastic to totally inelastic in 
four possible profiles. Given the values of the two measured 
strains, £ } and £ 2j the related profile is directly selected, and 
the zones and integration limits are decided as shown in 
figure 4. 



(I) e cy <£j<£2<£ry (II ) e C y<£l< E ty< £ 2 
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(Ill) ej <£ C y<£ty<£2 
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(IV) £ t y<£l<£2 


Fig. 4. Strain profiles over element cross 
section 

The induced combined normal and bending moment 
loads in each strain profile are resolved through a 
straightforward analytical routine summarized as follows: 

* Using azone boundaries strains from selected Fig. 
4 profile, the integration limits C aib for each zone are 
defined from Eq. (10), and are substituted into Eq. (9) to 
solve for the normal load of each zone in the 

profile, 



• The sum of normal loads from all zones, Eq. 
(11), is substiiuted into Eq. (5) to obtain the profile 
normal strain. 

• The bending neutral axis Cm is located using the 
total normal strain from Eq. (5) in Eq. (12) 

• Equation (14) determines the bending moment 
M lfi in each zone about the neutral bending axis, which 
are summed in Eq. (15) to provide the desired profile 
moment. 

• Strain distribution Ey and stress distribution a y 
are plotted over the thickness using Eqs. (6) and (16), 
respectively. 

This direct, though laborious, routine was reduced 
to a simple computer code requiring no detailed knowledge 
of its derivation. 

V. Normal-bending loads program 

Profile (III), having the most zones, was solved and 
programmed as summarized above. Other profiles, having 
fewer zones, were adapted by resetting limits according to 
their zone boundary values and positions in the strain 
diagrams and applying them to their appropriate zones. 

'NORMAL/BENDING LOADS FROM STRAIN DATA 
'NMLFSD, Microsoft Quick Basic 

' MATERIAL PROPERTIES 

INPUT "ELASTIC MODULUS E=";ELM 

INPUT "YIELD STRESS Fty=";FTY 

INPUT "MAX STRESS Ftu=";FTU 

INPUT "STRAIN @ MAX STRESS Etu=";ETU 

ETY=FTY/ELM 

PRINT "TENSION YIELD STRAIN";ETY 
ECY=-ETY 

SHE=LOG(FTU/FTY)/LOG(ETU/ETY) 

PRINT "STRAIN HARDENING EXPO. n=";SHE 

K=FTY/(ETY A SHE) 

PRINT "STRENGTH COEF K=";K 
KO=K 

SHE0=SHE 

ECY0=ECY 

ETY0=ETY 

'TEST DATA 

INPUT "ELEMENT THICKNESS H=";H 
INPUT "ELEMENT WIDTH w=";W 

10 INPUT "TEST MAX STRAIN E2=":E2 
INPUT "TEST MIN STRAIN E1=";E1 
IF E2<E1 THEN 

PRINT "MAX STRAIN < MIN STRAIN" 

GOTO 10 
END IF 

IF E2=E1 THEN E 1 =0.975 
SLOP=(E2-El)/H 


PRO=3 

USING PROFILE (HI) (E1<ECY<ETY<E2) 

IF ECY<E1 AND EI<ETY AND ETY<E2 THEN 
ECY=El:PRO=2 

ELSEIF ETY<E1 AND E1<E2 THEN 
EC Y=E 1 :ETY=E 1 :PRO=4 
ELSEIF E2<ETY AND ECY<E1 THEN 
K=ELM :SHE=l:ECY=El:ETY=E2:PRO=l 
END IF 

NIII1=W*K*(E2 A (SHE+1)- 
ETY A (S HE+ 1 ))/(SLOP” (S HE+ 1 )) 
NIII2=W*ELM*((ETY A 2)-(ECY A 2))/(2*SLOP) 
NIII3=(ABS(ECY)) A (SHE+ 1)-(AB S(E 1 )) A (S HE+ 1) 
NUI3=Nm3* W*K/(SLOP*(SHE+ 1)) 

niht=niii 1 +Nm2+Nm3 

PRINT "TOTAL AXIAL LOAD N=";NIIIT 
S NIII=NIIIT/W/H 

PRINT "AXIAL LOAD STRESS SN=";SNIII 

IF SNIIRFTY THEN 

ENIII=SNIII/ELM 

ELSE 

ENIH=(SNIII/K) A (1/SHE) 

END IF 

PRINT "AXIAL LOAD STRAIN EN="ENIII 
EMMIII=E2-ENIII 

PRINT "MAX BENDING STRAIN EM=";EMMIII 

CMin= (ENIII-El)/SLOP-H/2 

PRINT "BENDING NEUTRAL AXIS CM=";CMII 

Mil 1 =-((E2 A (SHE+ l))-(ETY A (SHEf 1 )))/(S HE+ 1 ) 
MIIl=MIIIl*((El+E2)/2+CMIIl , ‘SLOP) 

Mil 1=MII 1 +((E2 A (SHE+2))-(ETY A (SHE+2)))/(SHE+2) 
MII1=MIII1*W*K/(SL0P A 2) 

MII2=-((ETY A 2)-(E 1 A 2))*((El+E2)/2+CMlIRSLOP)/2 

MI12=MII2+((ETY A 3)-(El A 3))/3 

MII2=MI12*W*ELM/(SLOP A 2) 

MHI3=-(( AB S (EC Y)) A (SHE+ 1 )- 
( AB S (E 1 )) A (S HE+ 1 ))/(SHE+ 1 ) 
MII3=MIII3*((El+E2)/2+CMIII*SLOP) 
MII3=MIII3+((ABS(ECY)) A (SHE+2)- 
( AB S (E 1 )) A (S HE+2))/(S HE+2) 

MI 13 = M II 1 3 * W * K/(S LOP A 2) 

MIIT=MII 1 +MII2+MII3 

PRINT "BENDING MOMENT M=";M13T 

RIII=MIIT/NIIIT 

PR INT" MOMENT/AXI AL LOAD RATIO R=";RIII 
LIMITS 

CTY =(ETY -E l)/SLOP-H/2 
CCY=(ECY-El)/SLOP-H/2 
ETY A =FTY /ELM 

' STRESS & STRAIN DISTRIBUTIONS 
OPEN "CLIP:" FOR OUTPUT AS S2 

PRINT ”PROFILE=";PRO 
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IF PR0=3 THEN 

YS=-.5*H: YF=CCY : MY=9 

M=MY-1 

DY=(YF-YS)/M 

EY3=0: SY3=0 

y=YS 

FOR 1=1 TO M 

EY3=(.5*H+y)*SLOP+El 

SY3=K“((ABS(EY3) A SHE))*SGN(EY3) 

WRITE S2,y,EY3,ENIII,ETYA,SY3,SNIII,FTY 

PRINT yEY3,ENIII,ETYA,SY3,SNIII,FTY 

y=YS+(I+l)*DY 

NEXT I 

END IF 

IF PR 0=1 OR PR0=2 OR PR0=3 THEN 

YS=CCY: YF=CTY: MY=9 

IF E2<ETY THEN YF=.5*H 

M=MY-1 

DY=(YF-YS)/M 

EY2=0: SY2=0 

y=YS 

FOR 1=1 TO M 

EY2=(.5*H+y)*SLOP+El 

SY2=ELM*EY2 

WRITE #2,y,EY2,ENIII,ETYA,SY2,SNIII,FTY 

PRINT y,EY2,ENIII,ETYA,SY2,SNIII,FTY 

y=YS+(I+l)*DY 

NEXT I 

END IF 

IF PR0=2 OR PRO=3 OR PR0=4 THEN 

YS=CTY: YF=.5*H: MY=11 

M=MY-1 

DY=(YF-YS)/M 

EP1=0: SP1=0 

y=YS 

FOR 1=1 TO M 

EPl=(.5*H+y)*SL0P+El 

SPl=K*((ABiS(EPl) A SHE))*SGN(EPl) 

WRITE #2,y JEP1 ,ENIII JETY A,SP 1 ,SNIII FT Y 

PRINT y,EPl,ENIIIETYA,SPl,SNIILFTY 

y=YS+(I-rl)*DY 

NEXT I 

END IF 

CLOSE #2 
REM STOP 
CLS 

ETY=ETYO 

ECY=ECYO 

K=KO 

SHE=SHEO 
GOTO 10 


A sample printout of the program giving cross 
section characteristics derived from back-to-back strain gauge 
data is 

ELASTIC MODULUS E=? 10500 

YIELD STRESS Fty =? 38 

MAX STRESS Flu =? 58 

STRAIN @ MAX STRESS Eyu =? .06 

TENSION YIELD STRAIN 3.619048E - 03 

STRAIN HARDENING EXPO, n = .1505829 

STRENGTH COEF K= 88.59669 

ELEMENT THICKNESS H =? 1.4 

ELEMENT WIDTH w =? .74 

TEST MAX STRAIN E2 =? .02 

TEST MIN STRAIN El=? -.01 

TOTAL AXIAL LOAD N= 16.21604 

AXIAL LOAD STRESS SN= 15.65255 

AXIAL LOAD STRAIN EN= 1.490719E-03 

MAX BENDING STRAIN EM= 1.850928E-02 

BENDING NEUTRAL AXIS CM= -.1637665 

BENDING MOMENT M= 17.29161 

MOMENT / AXIAL LOAD RATIO R= 1.066328 

VT. Conclusions 

Experimental verification consists of two coherent, 
deterministic static test parts. Structural response within the 
elastic limit is verified with specified external loads 
representing maximum predicted operational environments. 
The ultimate factor of safely covers rare events, and its 
traditional and historical usage exerts the greatest influence 
on design and acceptance criteria. 

However, the ordcr-of-magnitude larger strains, and 
therefore displacements, imposed by the ultimate factor of 
safety may distort the applied load transmission. The 
documented technique was developed to idenufy and assess 
verification load transfer discrepancy through back-to-back 
surface mounted strain gauge data, which is applicable 
throughout the elastic and inelastic range of the structural 
material. 


It is concerning that verification test results often 
report surface strain measurements to conform very well 
with predicted math models up to the yield point, but then 
unexpectedly deviate during the inelastic loading to 
premature fracture. Reasons offered are usually indefinite. 
Perhaps this suggested technique may extend the basis for a 
more definite test evaluation. 
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